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ABSTRACT 



00 

^ Using single channel infrared photometry from the literature, a provisional K-band 

diameter-velocity dispersion relation for elliptical galaxies in the Coma and Virgo clus- 
\Q ters is derived. The Coma cluster relation has ^1.5 times lower scatter at K than in B 

CO or V. Excluding 4 outliers, the RMS scatter at K for 24 galaxies in Coma is only 4.8% 

in distance, close to the limit implied by the observational errors. Distance estimates 
based on the IR D — a relation will be more accurate than those derived from optical 
data. The improvement in the infrared is attributed to a decrease in sensitivity to 
stellar population parameters (age, metallicity, and slope of the IMF) as well as lower 
P { internal extinction from dust compared to the optical. That the D — a relation has a 

^ larger scatter in the optical indicates that there are detectable, but small, stellar popu- 

^ lation or dust content differences among the Coma ellipticals. Since the D — a relations 

are based on the fundamental plane, this result promises that the fundamental plane is 
thinner in the infrared than it is in the optical. 

Infrared photometric data available for Virgo are limited to just 13 objects; the 
spread in distance due to the depth of the cluster precludes any significant improvement 
over B and V. A relative Coma- Virgo distance of 5.56 is derived from the K band data, 
in agreement with estimates in other colors and using other techniques, indicating that 
there is no significant age difference between Virgo and Coma ellipticals. 



Subject headings: galaxies: distances — galaxies: elliptical 
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1. Introduction 

There are two popular methods for determin- 
ing redshift-independent distances to galaxies. 
Both rely on a tight correlation between a spec- 
troscopically determined distance-independent mea- 
sure of the gravitational potential and a photo- 
metrically determined distance-dependent mea- 
sure of apparent size. The Tully-Fisher (TF) re- 
lation for spirals (Tully & Fisher 1977) is based 
on the correlation between luminosity and rota- 
tional velocity. The D n — a technique for ellip- 
ticals (Dressier et al. 1987; D87) depends on the 
correlation between a surface brightness depen- 
dent diameter and the central velocity dispersion 
(a). 

When first introduced, the TF method em- 
ployed B band photometry, necessitating rela- 
tively large and uncertain corrections for inter- 
nal extinction which vary with galaxy morphol- 
ogy and star formation rate (Aaronson, Huchra, 
& Mould 1979 AHM). Differing current SFR's 
probably also contribute to the spread in B mag- 
nitudes. The accuracy of the TF method was 
improved markedly when AHM introduced the 
use of H band (1.6//) photometry where extinc- 
tion corrections are reduced to negligible levels 
and the flux is mainly from the oldest popula- 
tion of stars which is perhaps more uniform from 
galaxy to galaxy (AHM). 

There is increasing evidence from spectro- 
scopic studies suggesting that the stellar popu- 
lation of ordinary elliptical galaxies can be a mix 
of ages {e.g., O'Connell 1980; Rose 1985; Pickles 
1985, Bower et al. 1990, Gonzales 1993). Ex- 
treme cases such as merger remnants have expe- 
rienced strong bursts of star formation (Carter 
et al. 1988; Schweizer et al. 1990). Dust has 
also been shown to exist in a large fraction of 
ellipticals. There is good evidence that the dis- 
tances derived from D n — a and the closely related 
fundamental plane (FP) technique (Djorgovski & 
Davis 1987) for some ellipticals are influenced by 
one or more of these factors (Gregg 1992, 1995; 



Guzman & Lucey 1993). By analogy with the 
TF relation, it is natural to suspect that replac- 
ing B-band with infrared photometry might lead 
to an improvement in the accuracy of distances 
estimated by D n — a and FP plane techniques 
(Guzman & Lucey 1993; Gregg 1995). This pa- 
per presents evidence that the D — o method, and 
by inference the FP, has significantly less scatter 
in the IR than in the optical for ellipticals in the 
Coma cluster. 

2. Derivation of Infrared Diameters 

Using IR data from Bower, Lucey, & Ellis 
(1992a; BLE), it has been shown elsewhere that 
the K-band (2.2/^) magnitude- velocity dispersion 
relation for 25 early type galaxies in the Coma 
cluster has a scatter smaller by a factor of 1.9 
than the B band D — a relation (Gregg 1995). 
Because the IR magnitudes are from single chan- 
nel photometry in a 17" aperture, the derived 
IR relation is between velocity dispersion and a 
metric size. For a useful distance indicator, it 
is necessary to define a size based on structural 
parameters of the galaxies individually; the FP 
or D — a relations accomplish this using the half 
light-radius R e or D n . D87 define D n as the diam- 
eter of a circular aperture within which the mean 
B band surface brightness is 20.75 mag/sq arc- 
sec and show that log(D n ) is well correlated with 
log(cr), reducing the scatter by a factor of two 
from the Faber- Jackson (Bt vs. log(cr)) relation 
(Faber & Jackson 1976). Typical ellipticals have 
B-K ~ 4.0, so a K band diameter (Dk) which ap- 
proximates D n would enclose a mean K surface 
brightness of 16.75 mag/sq arcsec. 

The Virgo and Coma D n are based on CCD 
photometry (D87; Burstein et al. 1988, B88). 
The ideal IR version of a D — a relation would use 
K band imaging 1 ; however, with some simple as- 

1 Jorgensen, Franx, & Kjaergaard (1993) have demon- 
strated that the FP is an intrinsically more accurate dis- 
tance indicator than D — a. Using the fits detailed here, 
it is possible to derive an IR version of the FP. Because 
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sumptions a provisional Dk — o relation can be 
derived for the Coma and Virgo clusters from the 
data of Persson, Frogel, & Aaronson (PFA 1979) 
and BLE; both present high quality single chan- 
nel IR photometry for early type galaxies in the 
Coma and Virgo clusters. The BLE photometry 
has been calibrated to be on the same system as 
that of PFA (see BLE for discussion) and their 
much larger body of photometry supersedes the 
earlier PFA work. As BLE present no new IR 
photometry for Virgo, the data of PFA are used 
here. 

Table 5 of BLE lists K magnitudes through an 
effective aperture of 17" for 25 E or E/SO and 3 
SO galaxies in the Coma cluster which also have 
log(cr) > 2.0 in Davies et al. (1988) or Dressier 
(1987). Corrections for background light, red- 
shift, and Galactic reddening were applied as de- 
tailed in BLE and PFA. Dk can be derived from 
these single aperture data by assuming an R 1//4 
law profile and adopting an effective diameter 
A e from the V band photometry of Lucey et al. 
(1991a LGCT) for the ellipticals and the Third 
Reference Catalog (de Vaucouleurs et al. 1991) 
for the SO's. In the absence of color gradients, A e 
is independent of wavelength and in most ellipti- 
cals, color gradients are negligible (PFA; Peletier, 
Valentjin, & Jameson 1990). 

The R 1 / 4 fits are uniquely determined by the 
aperture photometry and effective diameter. Er- 
rors have been estimated by deriving Dk after 
varying both quantities by their estimated stan- 
dard deviations, ±0.028 dex for log{A e ) (LGCT) 
and ±0.026 mag for the K aperture magnitudes 
(BLE). These small errors have little effect on 
log(DK), amounting to only ±0.01 dex RMS, 



of the use single channel IR data and A e from V band 
measurements, this is not preferred here because the V 
band photometry then plays a role in determining 2 of 
the 3 dimensions of the FP (A e and /i e ), and contributes 
to both sides the FP distance estimator itself, A e as well 
as log (a ) + C*/x e . The results, however, of such an analysis 
are in agreement with the D — a relation discussed here. 
Several groups are now working towards a full K-band FP 
derivation (e.g., Pahre et al. 1995). 



smaller than the error estimates for log(D n ) Burstein 
et al. (1987) and comparable to those in log(Dy) 
(LGCT). The derived log(D K ) are listed in Table 
1 along with their A e , D n from B88, and the V 
band equivalent, Dy (aperture diameter enclos- 
ing 19.80 mag/sq arcsec) from BLE. Note that 
although V band A e from LGCT have been used 
in deriving Dk, strict photometric independence 
between the V and K band results is maintained 
by adopting Dy from BLE. 

PFA present JHK photometry of 13 Virgo 
ellipticals through 2 apertures (29" and 56"), 
which allows Dk to be derived from an R 1//4 law 
fit without recourse to any optical data. For 
comparison, Dk values were also derived for the 
Virgo ellipticals using the same method as the 
Coma ellipticals, i.e., by adopting an optical A e 
and fitting to a single aperture. It is testimony 
to the high quality of the data from PFA, BLE, 
and B88 that the two methods give identical re- 
sults within a few percent and that these are in 
turn extremely well correlated with the diame- 
ters in B and V. This test serves as a confidence 
check on the derived IR diameters for the Coma 
galaxies. Table 1 includes the 13 Virgo ellipti- 
cals from PFA and their calculated Dk, along 
with D n (B88) and Dy (BLE). One of the ob- 
jects, NGC 4387, with log(o-) < 2.0, is excluded 
from the analysis to follow. 

2.1. Velocity Dispersions 

Velocity dispersions from Davies et al. (1987) 
have been adopted because these data produce 
tighter D — a correlations than other data sets 
(e.g., Dressier 1984 or LGCT). This can be at- 
tributed to the careful treatment of systematic 
effects as well as the averaging of several inde- 
pendent measurements in many cases. In 3 cases, 
however, the velocity dispersions of LGCT are 
preferable. The LGCT a for RB 257 is more 
than 4 rms lower than that of Davies et al. and 
moves RB 257 directly onto the D — a relation 
in all colors. The Davies et al. dispersion for 
NGC 4874 is only 245 km/s, well below other es- 
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Table 1 

Data for Coma and Vireo Galaxies 



Name 


^ TT 


(j 


A e 


D n 


D v 


D K 






Coma 


Cluster 






NGC 4839 


31 


2.413 


1 79 ^ 

1. / zo 


1 9Q8 

1 .zoo 


1.278 


1.306 


NGC 4Q2fi 


49 


2.431 


1 981 


1 9/18 


1.269 


1.293 


TC 3959 

1 V J 'J <J fj <J 


69 


2.301 


l.UOU 


1 0/1 8 
1 .U4o 


1.064 


1.093 


TC 3957 


70 


2.176 


u.yzi 


u.yoo 


0.948 


0.932 


NGC 4923 


78 


2.318 


1 1/17 




1.144 


1.137 


NGC 4906 


118 


2.225 


1 1 /1Q 


i n^s 


1.029 


1.030 


NGC 4876 


124 


2.260 


u.yyu 


i nns 


1.030 


1.028 


RB 43 


125 


2.221 


n ^i ^ 


U.oDU 


0.874 


0.825 


NGC 4874 


129 


2.446 


9 0/10 
Z.U4U 


1 988 
1 .zoo 


1.287 


1.324 


NGC 4872 

_L > V 1 V V (J 1 _ 


130 


2.326 


0.768 


1.058 


1.048 


1.054 


NGC 4867 


133 


2.346 


0.852 


1.058 


1.060 


1.070 


RB 257 


136 


2.130 


0.446 


8(S8 
u.ouo 


0.877 


0.865 


IC 4051 


143 


2.348 


1.572 


1.078 


1.107 


1.134 


NGC 4889 


148 


2.581 


1.775 


1.468 


1.481 


1.544 


IC 4011 


150 


2.025 


0.976 


0.868 


0.868 


0.854 


NGC 4886 


151 


2.215 


1.173 


1.048 


1.036 


1.006 


IC 3998 


152 


2.192 


1.18 


0.895 


0.903 


0.934 


RB 45 


153 


2.130 


U.OIO 


U.OfJO 


0.842 


0.863 


IC 4045 


168 


2.324 


0.960 


1.098 


1.123 


1.148 


IC 4026 


170 


2.146 


1.16 


0.917 


0.921 


0.874 


IC 4021 


172 


2.199 


fl 777 

U. i i i 


Q^8 


0.956 


0.955 


IC 4012 


174 


2.253 


0.577 


0.978 


0.999 


1.032 


NGC 4883 


175 


2.234 


1.00 


1.026 


1.027 


1.050 


RB 155 


193 


2.079 


U.oOo 


818 
U.olo 


0.828 


0.808 


NGC 4860 


194 


2.391 


1.185 


1.168 


1.192 


1.228 


RB 167 


207 


2.167 


0.885 


0.898 


0.894 


0.894 


NGC 4881 


217 


2.340 


1.293 


1.108 


1.134 


1.146 


NGC 4841A 


240 


2.422 


1.459 


1.248 


1.299 


1.302 






Virgo 


Cluster 








NGC 4365 




2.394 




1.88 


1.88 


1.96 


NGC 4374 




2.458 




2.03 


2.01 


2.07 


NGC 4387 




1.924 




1.48 


1.47 


1.47 


NGC 4406 




2.398 




2.02 


1.94 


2.02 


NGC 4458 




2.025 




1.39 


1.39 


1.37 


NGC 4472 




2.458 




2.14 


2.14 


2.23 


NGC 4473 




2.250 




1.87 


1.86 


1.88 


NGC 4478 




2.174 




1.70 


1.67 


1.69 


NGC 4486 




2.558 




2.08 


2.09 


2.21 


NGC 4552 




2.417 




1.93 


1.90 


1.97 


NGC 4621 




2.381 




1.91 


1.90 


1.94 


NGC 4636 




2.281 




1.86 


1.87 


1.92 


NGC 4660 




2.296 




1.69 


1.71 


1.73 



timates (c.f. Whitmore, McElroy, & Tonry at 304 
km/s). The intermediate log(cr) of LGCT (279 
kms -1 ) again moves NGC 4874 onto the D — a 
relation in all bands. Finally, the Davies et al. 
value for NGC 4923 is the only one in the present 
sample based only on lower resolution Lick data; 
the higher precision LGCT value is superior. By 
adopting the alternative velocity dispersions for 
these three objects, the scatter in the D — a re- 
lation is minimized in all three colors. 

Additional support for these three velocity dis- 
persions being superior comes from the funda- 
mental plane: in planar fits to log(er) and /J, e of 
the other galaxies in B and V, the alternative dis- 
persions lower the residuals of the three objects 
in question, with the exception of RB 257 in the 
B band, by 2-5<r. The effective radius of RB 257 
varies by a factor of more than 1.5 in the liter- 
ature, ranging from 1.4" to 2.2". Seeing effects 
may be dominating the errors in R e for this ob- 
ject and its effective radius must be considered 
uncertain. 

3. Coma 

The IR and optical D — a relations for the 
Coma objects are compared in Figure 1 (upper 
panels, a-c). The filled circles are ellipticals and 
the 4-pointed stars are SO's. The open circles 
are 3 of the highest surface brightness galaxies in 
the sample (RB 43, NGC 4867, and NGC 4872). 
These clearly lie off the trend in K and tend to 
be high in V and B as well. They will be dis- 
cussed further below and are not used in the fits. 
The "+" is IC 4011 which has a log(er) just above 
the cutoff of 2.0 applied by D87; it lies well be- 
low the relation in all 3 colors and so will be 
excluded from further consideration. The lines 
through the data in Figure 1 are derived from 
a robust fitting technique based on median val- 
ues (Emerson & Hoaglin 1983), with log(cr) as 
the independent parameter, to the 24 remaining 
galaxies. Due to the nature of the robust fitting 
technique, there is little difference between fits 
with and without the outliers. The SO galaxies 
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Fig. 1. — Upper panels: comparison of the B, V, and K D — a relations for the Coma cluster. Filled 
circles are ellipticals, 4-pointed stars are SO's, open circles are outliers identified in the K band data (see 
text). The "+" is IC 4011. Lines are robust fits to the filled circles and 3 SO's. Lower panels: residuals 
of the 24 E and SO objects used in the fitting. The RMS in K is 1.68 times lower than in B, 1.42 times 
lower than in V. A representative error bar is plotted showing ±0.02 in Log(D). The error estimates for 
Log(DK) and Log(Dy) are half this amount. 



fit well among the ellipticals, as demonstrated 
by Dressier (1987). The best fits and Spearman 
correlation coefficients and probabilities are 

Log(D n ) = 1.279 log(cr) - 1.855; (1) 
S = 0.96; P s = 7.1 x 10" 14 

Log(Dy) = 1.319 log(cr) - 1.938; (2) 
S = 0.97; P s = 4.93 x 10~ 15 

and 

Log(D K ) = 1.458 log(cr) - 2.243; (3) 
S = 0.98; P s = 9.26 x 10 -18 . 



The residuals in Log(D) at a given velocity dis- 
persion are plotted in the lower panels of Fig- 
ure 1. 

Immediately apparent from Figure 1 is that 
the Dk — a relation is tighter than either of its 
optical counterparts. The RMS scatter of the 24 
galaxies is 0.0345, 0.029, and 0.0205 in log(D) for 
B, V, and K, respectively. The K band relation 
reduces the scatter by a factor of 1.68 from B 
and 1.42 from V. This corresponds to a scatter 
in distance of only 4.8%, roughly the predicted 
limit from measurement errors (Jorgensen et al. 
1993). The RMS for B and V are both less than 
that found by Jorgensen et al. from the Coma 
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FP for a similar sample; this can be attributed 
to the use of the alternative velocity dispersions 
and removal of the 3 remaining outliers. It can be 
demonstrated that removing the outliers and us- 
ing the better velocity dispersions also improves 
the FP fits. 

In K, there are three objects which lie well 
above the tight relation defined by the major- 
ity. At V, these objects are less clearly outliers 
and at B, only RB 43 is clearly away from the 
rest, due to the increase in sample RMS. The 
trend with bandpass suggests stellar population 
effects, but these three objects have typical V-K 
(BLE). These 3 are among the 5 highest surface 
brightness objects in the sample. Jorgensen et 
al. have shown that the correlation of residuals in 
the D — cr relation with surface brightness (Lucey, 
Bower, & Ellis 1991) results from the D — a rela- 
tion being a not-quite-optimum projection of the 
FP. While this may account for these outliers, it 
is noted that RB 257 and IC 4012 also have sur- 
face brightnesses in the same range, yet lie close 
to the D — cr relations in all colors (though not if 
log(cr) from Davies et al. is adopted for RB 257). 
LGCT redetermined log(cr) for RB 43 and it is 
in good agreement with the Davies et al. value. 
Unusual central dynamics may be responsible, 
but further longslit observations are needed to 
resolve this (Jorgensen et al. 1993). An indepen- 
dent, comprehensive spectroscopic study at high 
signal-to-noise and high dispersion is needed to 
place the Coma velocity dispersions on a truly 
uniform system. 

3.1. Statistical Significance 

The F-test indicates that the sample variance 
of Dk is significantly smaller than D n at the 
98.5% confidence level, while the improvement 
over Dy is significant at the 91.0% level. The 
precision of the photometry is such that the ve- 
locity dispersions are the limiting factor and con- 
tribute identically to the scatter in all three col- 
ors, resulting in highly correlated residuals. The 
F-test therefore underestimates the real signifi- 



cance of the differences. A Student's 2 point t- 
test takes this into account (Press et al. 1984) and 
indicates that the average absolute value of the 
residuals in Dk is significantly lower than in D n 
or D v at the 99.2% and 98.9% confidence levels, 
respectively. Either test applied between Dy and 
D n shows that they are not significantly different 
at the ~45% confidence level. The conclusion is 
that Dk — o is the best relation at a significantly 
higher level. 

3.2. Why is K better? 

Jorgensen et al. (1993) show that the observed 
scatter in the Coma FP in B and r is about twice 
the expected from the measurement errors, but 
could not identify a significant additional param- 
eter, ruling out ellipticity, B-r colors or gradients, 
isophotal twists, boxiness, and the presence of a 
disk. The photometric errors in B, V, and K are 
0.02-0.03 magnitudes; since the velocity disper- 
sions in Figure 1 are identical in B, V, and K, the 
increase in width of the D — cr relation with de- 
creasing wavelength is due to "cosmic" effects. 
One possible source of this additional scatter, 
which explains the decrease with wavelength, is 
variations of stellar population parameters (age, 
IMF slope, or metallicity) or the amount of dust 
(which may be linked to the stellar population) at 
a given velocity dispersion. All of these will have 
a larger effect in the optical than in the IR. Mod- 
eling carried out by Tonry, Ajhar, & Luppino 
(1990) shows that the I band is less affected by 
these variables than V or R, and it is reasonable 
to expect that the sensitivity is further reduced 
at K. Models by Bruzual (1983) suggest that this 
is the case (Guzman & Lucey 1993). The real 
situation may be a good deal more complicated, 
however. Early results from the surface bright- 
ness fluctuation distance indicator technique ap- 
plied in the K band show a larger scatter galaxy- 
to-galaxy than in the I band (Luppino & Tonry 
1993; Pahre & Mould 1994). If further observa- 
tions of nearby ellipticals continue to exhibit a 
large range in IR properties, the improvement in 
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the D — a relation demonstrated here for Coma 
in the K band may not be universal. 

The variations in age, IMF slope, and metal- 
licity can be quite small and still account for the 
observed increase in scatter. For instance, the 
total scatter in the B band D — a relation can 
be accounted for by a corresponding RMS in age 
of only 0.6-0.9 Gyr for simple models with mean 
age of 15 Gyr, typical IMF, and solar or twice 
solar metallicity (Gregg 1995, Table 1). Quite 
small changes in the metallicity or IMF can just 
as readily account for the scatter. This is an 
upper limit since some portion of the scatter is 
due to observational error rather than cosmic 
differences; taking the observational errors into 
account lowers the age scatter by 20-30%. The 
presence of dust would further reduce the allow- 
able scatter in the other stellar population pa- 
rameters. From these data alone, however, it 
is impossible to determine whether metallicity, 
age, the initial mass function, or simply dust is 
responsible for the additional scatter in the opti- 
cal. 

Similarly tight constraints on the ages have 
been reached by Bower, Lucey, & Ellis (1992b) 
using these same data in the form of color- 
magnitude relations. Although they detect a 
scatter in their Coma V-K color-magnitude re- 
lations larger than their estimated errors, they 
are reluctant to attribute it to a real variation in 
the Coma galaxies. The differences seen in the 
D — a relations of Figure 1 strongly argue that 
the larger scatter in the optical is due to intrinsic 
variations among the Coma early type galaxies at 
the tiny levels discussed above. 

Such tight constraints on the spread in the 
stellar population parameters stands in stark 
contrast to the often large range in ages at- 
tributed to ellipticals in the spectroscopic stud- 
ies cited in the Introduction. The disagreement 
can be resolved, as already suggested above, if 
the Coma cluster ellipticals contain a highly uni- 
form stellar population. Because most of the 
spectroscopic population studies have been con- 



ducted with narrow slits on relatively nearby ob- 
jects, aperture effects may also play a role. If 
the younger populations are confined to the in- 
ner 1-2" in galaxies 5-20 Mpc distant, they may 
be undetectable in Coma galaxies provided they 
have the same physical extent. A comparison of 
H/3 line strengths in the Coma and Virgo clus- 
ters (based on the same data as was used in the 
velocity dispersion determinations) by Gregg & 
Dressier (1986) shows that differences do exist 
and is consistent with an aperture effect. Bower 
et al. (1990), however, argue against aperture ef- 
fects as an explanation for the spectroscopic dif- 
ferences they find in a comparison of ellipticals 
from the Virgo and Coma clusters and other en- 
vironments. A much more detailed spectroscopic 
comparison of ellipticals in the Virgo and Coma 
cluster is called for to shed more light on this 
issue. 

Since the scatter in the D— a relation is smaller 
at K, then at a given luminosity or velocity dis- 
persion, the residuals in the optical D — a rela- 
tion should correlate with galaxy color. Galaxies 
which have too small a diameter in the optical 
at a given velocity dispersion must be brighter 
in the IR and hence redder for their total lumi- 
nosity; galaxies which are too large in the opti- 
cal will be bluer. Thus the diameter residuals 
from the D — a relation should anticorrelate with 
the residuals from the color-magnitude relation. 
Figure 2 plots the 24 residuals 5D Q of Figure 1 
against the residuals from the Vp — (V-K) rela- 
tion of Bower et al. (1992b), based on a sample 
of 40 E+S0 Coma galaxies. The Spearman rank 
correlation coefficient for the residuals is -0.63, 
probability 1.0 x 10~ 3 . The residuals in this plot 
are completely independent of each other and the 
evident anticorrelation strongly supports the re- 
sults of Figure 1 that the D — a relation is tighter 
in the IR and is consistent with the interpreta- 
tion that the additional scatter in the optical is 
from stellar population variations or dust. 

Jorgensen et al. (1993) report no correlation 
of residuals from a FP fit with B-r color resid- 
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the IR. Such work is already being carried out by 
several groups (e.g., Pahre et al. 1995). 
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Fig. 2. — Residuals 5D Q from the B band relation 
of Figure 1 plotted against residuals 5(V-K) from 
the V band color - magnitude relation of Bower 
et al. (1992b). The anti-correlation is expected if 
either stellar populations or dust are responsible 
for the increase in scatter of the B band D n — a 
relation. Symbols as in Figure 1. 

uals at a given velocity dispersion for the Coma 
galaxies. The positive result shown in Figure 2 
can be attributed primarily to greater sensitivity 
provided by the longer wavelength baseline of V- 
K, as well as complete independence of the two 
residual quantities. Correlated errors in residu- 
als both involving log(er) will tend to wash out 
any real effects from other causes. 

Even if the use of optical effective diameters to 
derive Dk for the Coma sample is inappropriate 
and the results of Figure 1 are completely fortu- 
itous, Figure 2 is an independent prediction that 
IR imaging data will show that the fundamental 
plane is thinner at K than in the optical. An 
IR imaging study of ellipticals in Coma, Virgo, 
and other clusters is needed for an independent 
assessment of the total improvement available at 
K and a complete characterization of the FP in 



4. Virgo 

Figure 3 compares the D n — a and D y — a rela- 
tions for 13 Virgo ellipticals with the Dk — o re- 
lation derived here, based on the PFA data. The 
lines are linear least squares fits, forcing the same 
slope as derived from the Coma relations, exclud- 
ing NGC 4387 ("+" sign) because of its low ve- 
locity dispersion. There is little difference among 
the colors and the RMS residuals are 2-3 times 
those for the Coma cluster. Within the uncer- 
tainties due to the large scatter and small num- 
ber of objects, the Virgo relations are compatible 
with those for Coma. The larger scatter of the 
Virgo ellipticals and the lack of improvement go- 
ing from B to K can be attributed primarily to 
the sizable depth of the Virgo cluster (e.g. Tonry 
et al. 1990). The computed intercepts for B, V, 
and K are -1.121, -1.224, and -1.498, respectively. 
From these, the relative Coma- Virgo distance in 
each color can be derived: 5.42 in B, 5.18 in V, 
and 5.56 in K, essentially identical given the scat- 
ter of ~0.07 in Log(D). These relative distances 
are in good agreement with other estimates (D87; 
Sandage & Tammann 1990; Bower et al. 1992b). 

Bower et al. (1990) present evidence for a 
rather large age difference between the ellipti- 
cal populations of the 2 clusters, based on Rose's 
(1985) spectral indices applied to composite spec- 
tra. The indicated age difference is approxi- 
mately a factor of two which in turn predicts 
that the derived distance ratio should be a rather 
strong function of wavelength, being smaller at 
longer wavelengths (Gregg 1995). No such trend 
is found, arguing for no significant difference in 
the mean ages of the Coma and Virgo ellipticals. 
A similar conclusion regarding the relative ages 
has been reached by Bower et al. (1992b) based 
on broad band colors of ellipticals in the 2 clus- 
ters. 
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Fig. 3. — Upper panels: comparison of the B, V, and K D — a relations for the Virgo cluster. The "+" 
is the low cr object NGC 4387. The lines are linear least squares fits (excluding NGC 4387), forcing the 
same slopes as in Figure 1 for the Coma cluster. Lower panels: residuals of the 12 E galaxies. The large 
scatter and high degree of similarity is attributed to the depth of the Virgo cluster. A representative 
error bar is plotted showing ±0.02 in Log(D). The error estimates for Log(Dx) and Log(Dy) are half this 
amount. 



5. Summary and Conclusions 

Using published single channel K band pho- 
tometry combined with optical effective diame- 
ters and the assumption of an R 1 / 4 law, a provi- 
sional IR version of the diameter — velocity dis- 
persion relation has been derived for early type 
galaxies in the Coma and Virgo clusters. The K 
band relation has 1.68 (1.42) times less scatter 
than the B (V) band D — a relation for 24 galax- 
ies in the Coma cluster, excluding a few outlying 
galaxies. The RMS scatter in distance is only 
4.8% per galaxy for the Coma cluster K band 
D — o relation, close to the limit from current ob- 



servational errors alone (Jorgensen et al. 1993). 

The improvement is attributed to lower sensi- 
tivity at K to variations in stellar population pa- 
rameters (age, metallicity, and slope of the IMF) 
as well as negligible influence from any dust. Dis- 
tance estimates based on K band photometry will 
be of superior accuracy to those in the optical. 
The improvement in the IR also indicates that 
there are small but detectable differences in ei- 
ther the stellar populations or the amount of dust 
in the Coma ellipticals. Because the D — a re- 
lation depends on the FP, a further implication 
is that the FP is thinner in the IR than in the 
optical. 
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A relative Coma- Virgo distance of 5.56 is de- 
rived from the K band relation, in good agree- 
ment with estimates from the B and V bands 
and other estimates from a variety of techniques. 
This argues against any significant systematic 
age difference between the elliptical galaxy pop- 
ulations of the two clusters. 

An IR imaging study of the Coma and Virgo 
ellipticals is needed to confirm the results derived 
here and to further investigate the FP at longer 
wavelengths. Additional velocity dispersion work 
at high resolution and S/N is needed to reduce 
the uncertainties in the velocity dispersion mea- 
surements which are the limiting factor in D — a 
and FP relations at all wavelengths. IR imag- 
ing and spectroscopy of other rich clusters would 
help to determine if the improvement seen in the 
Coma cluster K band D — a relation occurs in 
other environments and comparison with other 
colors could sensitively check the universality of 
D — cr and FP relations. 

This research has made use of the NASA/IPAC 
Extragalactic Database (NED) which is operated 
by the Jet Propulsion Laboratory, Caltech, un- 
der contract with the National Aeronautics and 
Space Administration. Matt Lehnert and Bob 
Becker are thanked for useful discussions and ad- 
vice. Part of this work was done at the Institute 
of Geophysics and Planetary Physics, under the 
auspices of the U.S. Department of Energy by 
Lawrence Livermore National Laboratory under 
contract No. W-7405-Eng-48. 
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